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Abstract 

In this paper a new lifetime distribution, which is called the exponentiated Weibull-geometric 
(EWG) distribution, is introduced. This new distribution obtained by compounding the ex- 
ponentiated Weibull and geometric distributions. The EWG distribution includes as special 
cases the generalized exponential-geometric (GEG), complementary Weibull-geometric (CWG), 
complementary exponential-geometric (CEG), exponentiated Rayleigh-geometric (ERG) and 
Rayleigh-geometric (RG) distributions. 

The hazard function of the EWG distribution can be decreasing, increasing, bathtub-shaped 
and unimodal among others. Several properties of the EWG distribution such as quantiles and 
moments, maximum likelihood estimation procedure via an EM-algorithm, Renyi and Shannon 
entropies, moments of order statistics, residual life function and probability weighted moments 
are studied in this paper. In the end, we give two applications with real data sets to show the 
flexibility of the new distribution. 

Keywords: EM-algorithm, Exponentiated Weibull distribution, Order statistics, Residual life 
function. 
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1. Introduction 

The Weibull and exponentiated Weibull (EW) distributions in spite of their simplicity in 
solving many problems in lifetime and reliability studies, do not provide a reasonable parametric 
fit to some practical applications. 

Recently, attempts have been made to define new families of probability distributions that 
extend well-known families of distributions and at the same time provide great flexibility in 
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modeling data in practice. One such class of distributions generated by compounding the 
well-known lifetime distributions such as exponential, Weibull, generalized exponential, expo- 
nentiated Weibull and etc with some discrete distributions such as binomial, geometric, zero- 
truncated Poisson, logarithmic and the power series distributions in general. The non-negative 
random variable Y denoting the lifetime of such a system is defined by Y = mini<j<jvXj or 
Y = maxi<j<jvXj , where the distribution of Xi belongs to one of the lifetime distributions and 
the random variable N can have some discrete distributions, mentioned above. 

This new class of distributions has been received considerable attention over the last years. 
The exponential-geometric (EG), exponential- Poisson (EP), exponential- logarithmic (EL), ex- 
ponential power series (EPS), Weibull-geometric (WG), Weibull-power series (WPS), exponen- 
tiated exponential Poisson (EEP), complementary exponential-geometric (CEG), two-parameter 
Poisson-exponential, generalized exponential power series (GEPS), exponentiated Weibull-Poisson 
(EWP) and generalized inverse Weibull-Poisson (GIWP) distributions were introduced and stud- 
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In this article, we propose a new four-parameter distribution, referred to as the EWG dis- 
tribution which contains as special sub-models the generalized exponential-geometric (GEG), 
complementary Weibull-geometric (CWG), complementary exponential-geometric (CEG), ex- 
ponentiated Rayleigh-geometric (ERG) and Rayleigh-geometric (RG) distributions. The hazard 
function of the EWG distribution can be decreasing, increasing, bathtub-shaped and unimodal. 
Several properties of the EWG distribution such as quantiles and moments, maximum likeli- 
hood estimation procedure via an EM-algorithm, Renyi and Shannon entropies, moments of 
order statistics, residual life function and probability weighted moments are studied in this 
paper. 

The paper is organized as follows. In Section 2, we review the EW distribution and its 
properties. In Section 3, we define the EWG distribution. The density, survival and hazard rate 
functions and some of their properties are given in this section. Section 4 provides a general 
expansion for the quantiles and moments of the EWG distribution. Its moment generating 
function is derived in this section. Renyi and Shannon entropies of the EWG distribution are 
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given in Section 5. Section 6 provides the moments of order statistics of the EWG distribution. 
Residual life function of the EWG distribution is discussed in Section 7. In Section 8 we explain 
probability weighted moments. Mean deviations from the mean and median are derived in 
Section 9. Section 10 is devoted to the Bonferroni and Lorenz curves of the EWG distribution. 
Estimation of the parameters by maximum likelihood via an EM-algorithm and inference for 
large sample are presented in Section 11. In Section 12, we studied some special sub-models of 
the EWG distribution. Applications to real data sets are given in Section 13 and conclusions 
are provided in Section 14. 



2. Exponentiated Weibull distribution: A brief review 



Mudholkar and Srivastava 
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which contains distributions with bathtub-shaped and unimodal failure rates besides a broader 



Gupta and Kundu [12|, Nassar and Eissa 



22 



class of monotone failure rates. One can see Mudholkar et al. |21], Mudholkar and Huston 

23j and Choudhury [9J for applications of the EW 



distribution in reliability and survival studies. 

The random variable X has an EW distribution if its cumulative distribution function (cdf) 
takes the form 

F x {x) = (l-e-^) a , x>0, (1) 

where 7 > 0, a > and (3 > 0, which is denoted by EW(a, f3, 7). The corresponding probability 
density function (pdf) is 

a-l 



fx{x) = a 1 p~<x~i- 1 e-^ x)1 (1 



-(£*)" 



(2) 



The survival and hazard rate functions of the EW distribution are 

S[x) = 1- [l-e-^y , 



and 



h{x) = a 7 /3T^- 1 e-^ 7 f 1 



a-l 



l_(l_ e -0*r 



respectively. The kth moment about zero of the EW distribution is given by 



E(X k ) = a/3- k T[- + 1 ) ]T(-1 
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Note that for positive integer values of a, the index j in previous sum stops at a — 1, and the 
above expression takes the closed form 

E(X k ) = ar k r(- + l]A k ( 1 ), (4) 
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where 

a-1 



A k ^) = l + J2(-iy(^ . ^( i + l)-(f +1 ), k = 1,2,3,-.. , (5) 
in which T(.) denotes the gamma function (see, Nassar and Eissa [23] for more details). 



3. The EWG distribution 

Consider the random variable X having the EW distribution where its cdf and pdf are given 
in © and ©. 

Given N, let Xi,--- ,Xjy be independent and identically distributed (iid) random variables 
from EW distribution. Let the random variable N is distributed according to the geometric 
distribution with pdf 

P(N = n) = (1 - 9)9 n ~ 1 , n = 1, 2, • • • , < 9 < 1. 

Let Y = max(Xi, • • • ,Xn), then the conditional cdf of Y\N = n is given by 

FY\N=n(y) = {l-e-^) na , (6) 

which is the EW distribution with parameters na, (3, 7, and denoted by EW(na, f3, 7). The 
exponentiated Weibull-geometric (EWG) distribution, denoted by EWG(a, f3, 7, 9), is defined 
by the marginal cdf of Y, i.e., 

Fy(y) = i-g(i- e -(^r • ( } 

This new distribution includes some sub-models such as the complementary exponential-geometric 
(CEG), generalized exponential-geometric (GEG), complementary Weibull-geometric (CWG), 
exponentiated Rayleigh-geometric (ERG) and Rayleigh-geometric (RG) as special cases. The 
pdf of the EWG distribution is given by 

y ) = r ; 72 ^v ai 2 ' 8 

[l-9(l-e-^) ] 
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where a, (3, 7 > and < 8 < 1. 

The survival function and hazard rate function of the EWG distribution, are given respectively 
by 

i_ (i- e -{pyyy 

S ^ = i_0(i_ e -(/3^) a ' (9) 

and 

(1 - B)ar l Fyt- 1 e-VvV (l - e-W) " 1 



[l-Q(l- e-^Y] [1 - (1 - e-(^) 7 )"] ' 

Proposition 1. For \9\ < 1, it is easy to prove that the density of EWG distribution can be 
written as an infinite mixture of EW distributions. If \z\ < 1 and k > 0, we have the series 
representation 

3=0 y u 

r 1 — 2 

If |0| < 1, expanding 1 — 6> (l — e"^) 7 )" as in Eq. (|10p . the density function ([5]) can be 
demonstrated by 

j + 1)V J II - e~^ y > 

.1 

Using the EW density ([2]), we obtain 



00 

f Y ( y ) = (l _ ^c^V^e-^ 7 j^O' + (* - e" (/3s/)7 )' 



/wc(!/;a,M,«) = (l-^^/w(!/;tt(i + l),/3,7). (11) 

3=0 

Various mathematical properties of the EWG distribution for \9\ < 1, can be obtained from Eq. 
(jlip and the corresponding properties of the EW distribution. 



Proposition 2. The density of EWG distribution can be expressed as infinite linear combi- 
nation of density of the biggest order statistic of X\,--- ,X n , where Xi ~ EW(a, /S, 7) for 
i = 1, 2, • • ■ , n. we have 



/Wy) = ^2(G(y)) n P(N = n) = £ 5X(n) (y)P(iV 
n=l n=l 

in which (y) is the pdf of X (n) = max(Xi, • • • , X n ). 
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4. Quantiles and moments of the EWG distribution 



The pth quantile of the EWG distribution is given by 

x P = /r 1 



log 1 



i-e(i-p) 



which is used for data generation from the EWG distribution. In particular, the median of the 
EWG distribution is given by 



xo.5 = & 1 



log 1 



1 



2-e 



Suppose that Y ~ EWG(a, f3, 7, 9), and Xi n \ = max(Xi, • • • , X n ), where Xj ~ EW(a, /3,"f) 
for i = 1,2, • • • , n, then the A;th moment of Y is given by 

E(Y k ) = E(E(Y k \N)) = EZi P(N = n)E(Y k \n) = ^Zi P(* = n)E(X k n) ) 

= EZi P(N = n)n«r fc r (A + l) E^o(-l) i CT') U + 

= (1 - 9)ap~ k r (4 + 1) xr=i ET=o^ n ~H-^ j (T 1 )^' + 

(12) 

For positive integer values of a, the index j in above expression stops at a — 1. 
Using Eq. f)12|) , the moment generating function of the EWG distribution is given by 

MY(t) = ZZo t i E ( Yi ) 

= E,=o I [(i - *)a/rT + 1) E~=i E^o^'C-iy CT 1 )^" + 1 )~ ( * +1> " 



a(i - o) ESo E£i Er= nfl-H-iy (T'K ( 7 + 1 



(i + 1) 



(13) 



According to Eq. (|12p . the mean and variance of the EWG distribution are given respectively 
by 



E(Y) = (1 - B)apr X T Q + l) E E ^("^ (™ ') 



and 



(„ s OO OO 



f ntt - l )(i + i)-- 



£ 2 (y), (15) 



Where is given in Eq. ([I 
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5. Renyi and Shannon entropies 

For a random variable with pdf /, the Renyi entropy is defined by Ir(t) = log{/ K f r (y)dy}, 
for r > and r / 1. For the EWG distribution, the power series expansion gives 

Jo°° f r (v)dy = [a7/? 7 (l - 0)T ET=o Jo°° yb-^e-rW (1 - e -(W)«CJ-H-)-r 

= [07^(1 -*)] r £^x^(-i)W 

Setting u = (/3y) 7 , gives 

Jo°° f r (y)dy = Hi - o)] r {fiif- 1 ET=o XXaC-^Y^ r^? T(r 'SL 0j ■ ( 16 ) 

Substituting from (116j) . the Renyi entropy is given by 

(17) 

The Shannon entropy which is defined by E[— \og( f(Y))], is derived from lim,,^! Ir(t). 

6. Moments of order statistics 

Let the random variable Y r:n denotes the rth order statistic (Y\ :n < Y<2- n < • • • < Y n - n ) in a 
sample of size n from the EWG distribution. The pdf of Y r:n for r = 1, • • • , n, is given by 

= ^ — — -/(y)F(y)'- 1 [1 - , y > 0. (18) 

-£>(r, n — r + 1) 

where F(y) and f(y) are the cdf and pdf of the random variable Y. Substituting from ([7J and 
© into (USD, gives 

f f./| - <*7/H(l-flr ■■ 7 -l--C8y)T (l-e-(^) 7 )— i[i-(i- e - W>)°]"- r 

Also the cdf of Y r:n is given by 

Fr:n{y) = ELr (*) TO)]* t 1 ~ W" * 

n (l- e) fe (l-e^^)" fc (l-(l- e -^ 7 ) a )"- fc 
- Z^fc=r (l-0(l-e-(/3f) 7 ) Q )" • ^ ^ 

Using the binomial expansion, series expansion (|10p and after some calculations, the fcth 
moment of the rth order statistic Y r:n is given by 

ap- k (l-d) r ^ j+sol (n + i\(n-r\(a(i+j + r)-l\ r(& + 1) 

l 



^)-c^5§5<-^( n n( B j r )r +i r ,_i ), +l) , 

(21 
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The pdf, cdf and kth. moment of the smallest and biggest order statistics, i.e., Y\ :n and Y n:n , can 
be obtained by setting r = 1 and n in Eqs. (I19| )- (l2ip . 

7. Residual life function of the EWG distribution 

Given that a component survives up to time t > 0, the residual life is the period beyond t 
until the time of failure and defined by expectation of the conditional random variable X\X > t. 
In reliability, it is well known that the mean residual life function and ratio of two consec- 
utive moments of residual life, determine the distribution uniquely (Gupta and Gupta, [111]). 
Therefore, we obtain the rth order moment of the residual life via the general formula 



m r (t) = E[(Y- t) r \Y > t ] = * jHfo - t) r f(y)dy, 



(22) 



where S(t) = 1 — F(t), is the survival function. 
In what seen this onwards, we use the expressions 



oo 

X 



J+s -l e - ik+ W x F dx = -JL_(jfc + ir (i + i) r(fc+1 )(^ (1 + £ }) 



and 



where r*(s) = x s ~ 1 e~ x dx is the upper incomplete gamma function and I^(s) = x s ~ 1 e~ x dx 
is the lower incomplete gamma function. 

Applying series expansion (flQ|) . the binomial expansion for (y — t) r and substituting S(y) given 
by @ into (|22p . the rth moment of the residual life of the EWG is given by 

(23) 

Another important representation for the EWG is the mean Residual life (MRL) function obtain 
by setting r = 1 in Eq. (|23p . The importance of the MRL function is due to its uniquely 
determination of the lifetime distribution as well as the failure rate (FR) function. Lifetimes 
can exhibit IMRL (increasing MRL) or DMRL (decreasing MRL). MRL functions that first 
decreases (increases) and then increases (decreases) are usually called bathtub (upside-down 
bathtub) shaped, BMRL (UMRL). Many authors such as Ghitany [loj], Mi 18], Park 24] and 



Tang et al. [26f] have been studied the relationship between the behaviors of the MLR and FR 
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functions of a distribution. 

The following theorem gives the MRL function of the EWG distribution. 
Theorem 1. The MRL function of the EWG distribution is given by 



mi(t) 



PS(t) 



YYi-iy ( a{l + + i)(j + i)-( i +^) r c?+ i )W(i + 1) 

V i / 7 



Setting r = 2 in (|23p . the variance of residual life function of the EWG distribution can be 
obtained using mi(t) and rri2(t). 

8. Estimation and inference 

In this section, we study the estimation of the parameters of the EWG distribution. Let 
Y±, Y2, • • • ,Y n be a random sample with observed values yi,y2, • • • , 2M from EWG distribution 
with parameters a, /3, 7 and 0. Let = (a, f3,~/,6) T be the parameter vector. The total log- 
likelihood function is given by 

In = ln(y, 8) = n[log a + log 7 + 7 log (3 + log(l - 0)] + (7 - 1) £™ =1 log y, - £? =1 C%) 7 

+ («-!) EJLi l°gU " e" (fe)7 ) - 2E?=i l«g[l " * (1 " e-(^) T ) a ]. 

The associated score function is given by U n {Q) = (dl n /da, dl n /d/3, dln/dj, dl n /d\) T , where 

= 5 + Er=i mi - e-^) 7 ) + 2^ Er=i io g (i - e -w p-*-^ 7 ?. 



l-0(l-e-(' 9 !' l ) 7 ) Qi ' 

§ =f -7^- 1 Er= 1 y7 + («-i)7^- 1 Er=iSS^ 



+ ZC, «7P 2vi=l i_fl(i_ e -(i8»i)'i , )o 



«* = 2 + n iog/3 + E"=i logy,. - EILi log(/3 yi )(^) 7 

, f„ ,\V» log^K^i)]^-^^ , ofl„,y^n log(/3 ; ; I )(/3^)^ e -(^) 7 (l- e -(' 9 ^) 7 )"- 1 

+ l« - 1 J 2^=1 1 _ e -(/3 3 /iJT r 2t/a Z^i=l i_ e( i_ e -(/3 ! / i )T) a ' 

9^ n_, 9 V- n (l-e-(^i) 7 )° 

~~ 1-0 "t z 2^i=i i_e(i_ e -(/3Hi)' y )a • 

The maximum likelihood estimation (MLE) of 0, say 0, is obtained by solving the nonlinear 
system U n (0) = 0. The solution of this nonlinear system of equation has not a closed form. For 
interval estimation and hypothesis tests on the model parameters, we require the information 
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matrix. The 4x4 observed information matrix is 



Icea 


la/3 




IaO 


IaP 






1/39 






J 77 




IaO 


1/36 




lee 



whose elements are given in Appendix. 

Applying the usual large sample approximation, MLE of i.e. can be treated as being 
approximately A^G, J n (0)~~ ), where J n (0) = E[I n (Q)]. Under conditions that are fulfilled 
for parameters in the interior of the parameter space but not on the boundary, the asymptotic 
distribution of ^/n(G-6) is A 4 (0, J(e) -1 ), where J (G) = lim 

n— >-oo n -^n(®) is the unit in- 
formation matrix. This asymptotic behavior remains valid if J(Q) is replaced by the average 
sample information matrix evaluated at G, say n^ 1 I n (Q). The estimated asymptotic multivari- 
ate normal A4 (0,7^(0) ) distribution of G can be used to construct approximate confidence 
intervals for the parameters and for the hazard rate and survival functions. An 100(1 — 7) 
asymptotic confidence interval for each parameter Q r is given by 

ACI r = (G r - Zq_ \fW\ Q r + Z2 \TF r ), 

where I rr is the (r, r) diagonal element of I n (©) for r = 1, 2, 3, 4, and Zx is the quantile 
1 — 7/2 of the standard normal distribution. 



9. Submodels of the EWG distribution 

The EWG distribution contains some sub-models for the special values of a, (3 and 7. Some 
of these distributions are discussed here in details. 

9.1. Complementary Weibull- geometric distribution 

The complementary Weibull-geometric (CWG) distribution is a special case of the EWG 
distribution for a = 1. Our approach here is complementary to that of Barreto-Souza et al. jj] 
in introducing the Weibull-geometric (WG) distribution. The pdf, cdf and hazard rate function 
of the CWG distribution, are given respectively by 

= K± "HP V e (24) 

[I - 9 {1 - e-(M")] 2 
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, , (l-«)(l-e-<»") , , 

F ^ = !-»(!-«-<*>■) ■ <25> 

and 

(1 - ^)7/3 7 y 7_1 e _(/3?/)7 

{y} ~ [i - e (i - e-wvF)] [i - (i - e-^yy)] ' (26) 

According to Eqs. (fTl)) and (fT5|) . the mean and variance of the CWG distribution are given by 

E(Y) = (l + I) fjEnr-H-lV 7 *) (J + (27) 



and 

(1 



Var(Y) 



r(i + ^EE nr-^-iy ( n " x ) (j + i)~^ - £ 2 (n (28) 



£2 



where -E(Y) is given in Eq. (1271) . 

9.5. Generalized exponential-geometric distribution 

The generalized exponential-geometric (GEG) distribution is a special case of the EWG 
distribution for 7=1. This distribution is introduced and analyzed in details by Mahmoudi and 
Jafari [15]. The pdf, cdf and hazard rate function of the GEG distribution, are given respectively 
by 



f(y) 



(1 - 9)a(3e-^ (l 



e 



[1 _ (1 _ e -m) a ] 2 

(l- g )(l-e-C»>)° 
W i-9(l-e-(M) a ' 

and 

(l-^e-^ri-e-^)) 01 " 1 

h(y) ~ 



[1-0(1- e ~(M) a ] [1 - (1 - e-(M) a ] ' 
According to Eqs. (|14p and (|15p . the mean and variance of the GEG distribution are 

n=l j=0 ^ 



and 

00 00 



Var(Y) = ?2^fi ^^^-(-ly ("7 l ) (J + I)" 3 - EHY), 



n=l j=0 

where E(Y) is given in Eq. (|29p . 
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9.3. Complementary exponential- geometric distribution 

The complementary exponential-geometric (CEG) distribution is a special case of the EWG 
distribution for a = 1 and 7 = 1. The pdf, cdf and hazard rate function of the CEG distribution 
are given respectively by 

(1 - 0)Pe-(M 
I[y) [l-0(l-e-^))] 2 ' 

(1 - 6) (1 - e-^) 
1 _ (i _ e -(M) ' 

and 

(1 - 6)(3e-^ 



F{y) 

Kv) 



[l-0{l-e-W)] [l-(l-e-W))]' 
According to Eqs. (|14p and (|15p . the mean and variance of the CEG distribution, are given 

respectively by 

BCn = (i-fl)/r 1 EEn^ 1 (-iy( n T /(j + i) 2 , (so) 

71=1 j=0 \ J / 

and 

OO OO / 1 \ 

Var(Y) = (1 - #)/T 2 r(3) E E (~ )/U + I) 3 - 

71=1 J=0 \ J J 

where -E'(^) is given in Eq. (I30p . 

9.^. Exponentiated Rayleigh- geometric distribution 

The exponentiated Rayleigh-geometric distribution (ERG) distribution is a special case of 
the EWG distribution for 7 = 2. The pdf, cdf and hazard rate function of the ERG distribution 
are given respectively by 

2(1 - 6)aP 2 ye-^ 2 (l - e^T" 1 

My) = — — „ NCn 2 . 



Fy(y) 



[l_^(l_ e -(/3y) 2 ) Q ] 
(1 - 0) ( 1 - e-Cw) 



l_0(l_ e -(^) 2 ) Q ' 
and 

2(1 - 0)a/3V- (/3y)2 (l-e-^) 2 )' 



%) 



[1 - (1 - e -(/3y) 2 ) a ] [1 - (1 - e -(W 2 ) a ] ' 
The mean and variance of the ERG distribution, are given respectively by 

e(y) = (i - e) a pr x v Q + i)EE ^-H-i)^ ( na ~ x ) (j + 1 



(1+1) 



(31) 



:r- 
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and 



Var(Y) = (1 - e)af3- 2 ^^^(-iW™ 7 ) (j + 1)~ 2 - E 2 (Y), 



n=l j=0 

where E(Y) is given in Eq. (I31|) . 

9.5. Rayleigh- geometric distribution 

The Rayleigh-geometric distribution (RG) distribution is a special case of the CWG distri- 
bution, obtained by choosing 7 = 2 in CWG distribution. Setting 7 = 2 in Eqs. (I24|) -(|26p gives 
the pdf, cdf and hazard rate function of the RG distribution as 

2(1 - e)l3 2 ye-^ 2 



f(y)-- 
F(y) 



[l_#(l_ e -W)] 2: 

(1 - e) (1 - e-(^) 2 ) 



l-6(l-e-(M 2 ) 

and 

2(1 - e)(3 2 ye 



2 \ ! 



2 7 , P -(/3y) 2 



e -{PyY [1 _ (1 _ e -(/^) 2 )] ' 
According to Eqs. (|27p and (|28p . the mean and variance of the RG distribution are given by 



9+ 1 EE" ri (- i ) j r • (32) 

Z ' n=l j=0 \ 3 J 



and 

00 00 



Far(F) = (1 - 0)/3" 2 £ £ ^ n_1 (-l) j ( " ~ )/U+ I)" 2 " E 2 (Y), 

n=l j=0 \ 3 J 

where E{Y) is given in Eq. ()32p . 
10. Conclusion 

We propose a new four-parameter distribution, referred to as the EWG distribution which 
contains as special sub-models the generalized exponential-geometric (GEG), complementary 
Weibull-geometric (CWG), complementary exponential-geometric (CEG), exponentiated Rayleigh- 
geometric (ERG) and Rayleigh-geometric (RG) distributions. The hazard function of the EWG 
distribution can be decreasing, increasing, bathtub-shaped and unimodal. Several properties of 
the EWG distribution such as quantiles and moments, maximum likelihood estimation proce- 
dure via an EM-algorithm, Renyi and Shannon entropies, moments of order statistics, residual 
life function and probability weighted moments are studied. Finally, we fitted EWG model to 
two real data sets to show the potential of the new proposed distribution. 
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